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Abstract 

Bogomolny's formula for energy-smoothed scars is applied for the first time 

to a non-specific, non-scalable Hamiltonian, a two-dimensional anharmonic 

oscillator. The semiclassical theory reproduces well the exact quantal results 

over a large spatial and energy range. 
PACS numbers: 03.65.Sq, 05.45.+b 
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The scars left on stationary state wavefunctions by unstable classical periodic orbits were 
discovered by McDonald and Kaufman [1J] . They were named and further studied by Heller 
. The paradox of the scars is that they occur in an energy region where a generic classical 
trajectory covers the energy shell uniformly; in spite of this, many quantal wavefunctions like 
to concentrate in the vicinity of one or several periodic orbits. Obviously, an understanding 
of this paradox requires a semi-classical theory. 

Such a semi-classical theory of scars in the usual, coordinate space wave functions was 
given by Bogomolny [[|. There are also semi-classical theories of scars in the Wigner distri- 
bution 0] and in the Husimi distribution ||. We present here some results from an extensive 
comparison of the exact quantal scars with the Bogomolny scars || for an ordinary smooth 
potential. Apart from an interesting qualitative discussion |J , we do not know of any other 
detailed comparison for a sufficiently "generic" Hamiltonian. 

The famous Gutzwiller trace formula [0] can be obtained by integrating the Bogomolny 
formula over the space coordinates. Therefore our work leads also to a check of the Gutzwiller 
trace formula for our Hamiltonian. There have already been a few checks of similar quality 
for the Gutzwiller trace formula with general enough Hamiltonians, for instance Ref. ||. 

Our Hamiltonian is 

H(p x , p y , x, y) = \(pl+ P 2 y ) + 0.05x 2 + (y - x 2 /2) 2 . (1) 

Its classical dynamics and periodic orbits have been studied in detail ||. We had to continue 
the study of periodic orbits towards higher energy, using a totally new method for which 
there is no space here. Suffice it to say that this method ensures that no orbit below period 
20 is missed. Fig. |l] shows a contour plot of the potential V(x, y) = H — |p 2 for energies of 
interest. The transition from mostly regular to mostly chaotic motion happens for energies 
of order 0.1. Islands of regularity remain, however, no matter how high in energy one goes. 
But at the energies of our numerical comparisons, they are very tiny. 

We took h = 0.05. The calculation of the exact quantal wave functions and energies 
was done with the basis cj)!^ ) (x)cj)^(y ~ x2 /2), where tfiffi and <p£i are harmonic oscillator 
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wave functions appropriate for the bottom of the well. We used 240 oscillator states for 
the x direction and 26 for y. A sensitive test of the sufficiency of this basis comes from a 



comparison with the smooth Thomas-Fermi density of states, including the corrections (TC 
of order h 2 . The test shows that our basis begins to fail for an energy somewhere between 
0.80 and 0.85, and consequently we stopped our comparisons there. 
Bogomolny's formula is || 

A(q,£; e )=— * £ fjr) ' ! 



(2^) 3 / 2 pc ^ d J eK \qi\ J\ 
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A is the oscillating part (as explained in the next sentence) of the energy-smoothed coordi- 
nate space probability density 

J2fe(E-E n ) |<q|^>| 2 , (3) 

n 

where f € is the smoothing function, which we take to be the normalized gaussian 

f e (E) = {2^h-\- E2 ^ 2 . (4) 

A is obtained by subtracting from Eq. the energy-smoothed Thomas-Fermi density. On 
the right hand side of Bogomolny's formula, f e (r) = e- V2 /2ft 2 is the Fourier transform of 
f t {E), r is the period of the periodic orbit, and S is its action / p • dq. The coordinates q\ 
and q 2 are chosen especially for each periodic orbit, qi being the distance along the orbit 
and q 2 being the perpendicular coordinate. M. is the 2x2 submatrix of the monodromy 
matrix involving coordinates q2 and P2, and m qp is one of its off-diagonal elements. Finally 



/i is equal to fx m , the Maslov index of the orbit fliT] , when m qp and TrA^ — 2 have the same 
sign; /i is equal to \i m — 1 when the signs are opposite. 

In order that Eq. (Q) be valid, it is also necessary to perform some smoothing over 
coordinate space on both sides of the equation. There are two reasons for this: (a) it is 
essential for reducing the semiclassical contribution to a sum over periodic orbits ; (b) the 
semiclassical theory is not valid near the points where m qp = 0, which are the self-conjugate 



points, and the spatial smoothing minimizes this discrepancy. We smoothed with a gaussian 
proportional to e~^ +q y^ b2 . This is equivalent to calculating the coordinate space projection 
of the Husimi distribution. We chose b = 0.2; more about this later. 

Gutzwiller's trace formula 0, obtained by integrating Eq. (|2]) over space, is 

d ( E ' £ ) = ^ J2 fe( T ) r—^ = C0S (f - (5) 

^pcoodic J\trM-2\ n 2 

orbits V 

d is the oscillating part of the energy-smoothed density of states, calculated by subtracting 
the Thomas-Fermi density of states and its corrections |10| of order ft 2 . For orbits consisting 
of repeated traversals of a primitive periodic orbit, r Q is the period of the latter. 

It has long been everybody's dream to get rid of the energy smoothing; to let e — > and 
to use Eqs. @ and (§) to predict individual stationary states. Like everyone else we avoided 
this limit, as the process seems to diverge. Instead, we chose our e so that relatively few 
periodic orbits would contribute. 

Fig. shows pictures of the first 12 orbits in order of increasing period, again at E — 0.8. 
We have worked in the range 0.5 < E < 0.85 and , for E ^ 0.8, the order might be slightly 
different from that in Fig. 0. In the distribution of periods, the three lowest (4.44, 6.44, 7.14) 
are clearly separated from those above (10.51, 11.57, 11.60, . . .etc.). Hence we expect that, 
by choosing e large enough, many features can be described in terms of 3 periodic orbits 
only. This turns out to be true indeed. Fig. |3] shows the oscillations in the smoothed energy 
level density, as given by the Gutzwiller formula, Eq. ([5]), calculated with e = 0.01 and 5 
periodic orbits, but it is only very slightly better than that calculated with 3 orbits. On the 
other hand it is radically different from that calculated with only 1 or 2 orbits. In the case 
of 3 or more orbits, the agreement with the exact quantal calculation can be termed very 
good for this value of e. 

Fig. |] shows the most striking scar we have found: odd (in x) state no. 145, with E = 
0.814, scarred by the third periodic orbit, a simple asymmetric libration. Our calculation 
of energy-smoothed scars according to the Bogomolny formula, Eq. (^), is illustrated in 
Figs. [|, ^, and [7] and compared with the exact result. For lack of space, we concentrate 



here on the energy range 0.726 < E < 0.814, which is typical of other energies. Fig. |5| 
is for the same energy as Fig. ^, and the scar is still prominent. Figs. |6| and [7] are two 
other examples. They all show some combination of scarring by the first three orbits. The 
agreement between the exact quantal and the semi-classical pictures can be described as 
following the main trends very well, but quantitatively inaccurate. By "main trends", we 
mean in particular the way that the density fluctuation in the vicinity of each orbit (the 
"scar" of that orbit) oscillates and changes sign as a function of energy. This is better seen in 
Fig. D, which shows the scar strength (both exact and semiclassical) as a function of energy 
for three points in the xy plane. 

The agreement between Eq. ([|) and the exact quantal density is limited to the central 
part of our potential, exclusive of the two "arms" (see Fig. |]). This is because, in the arms, 
the motion is approximately integrable, consisting of fast transverse oscillations whose action 



is an adiabatic invariant for the slower longitudinal motion fl2| . As a result, the exact density 
in the arms exhibits a very simple pattern, visible in the upper corners of Figs. |^, [(], and [7], 
and consisting mostly of wide longitudinal oscillations whose wavelength can be calculated 
very simply with the above adiabatic picture. The short-period classical orbits, on the other 
hand, are found only in the central, chaotic region. It is possible to increase the value of b, 
the spatial smoothing width, so as to wash out the oscillations in the arms, but then this 
would also wash out many interesting features in the center. The value chosen, b = 0.2, is 
a reasonable compromise. 
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FIGURES 



FIG. 1. Three contours of the potential V(x,y): E = 0.1 (dots), E = 0.4 (dashes), E = 0.8 (full). 

FIG. 2. The 12 orbits with lowest periods for E = 0.8. 

FIG. 3. The smoothed energy level density, with e = 0.01, minus the Thomas-Fermi den- 
sity. The full line is the exact quantal calculation, the dashed line is Eq. (5) with 5 orbits. The 
Thomas-Fermi density at E = 0.65 is ~ 290. Hence the average number of stationary states inside 
one e is 2.9. 

FIG. 4. Contours of the exact coordinate space probability density for odd state no. 145, 
showing striking scarring by the no. 3 periodic orbit and its mirror reflection x — > —x. 

FIG. 5. The oscillating part , A, of the energy-and-space-smoothed coordinate space density. 
Energy smoothing parameter e = 0.01. Top: exact quantal. Bottom: Eq. (2) with 5 orbits. Solid 
contours are positive, dashed are negative. The contour spacing is 0.5 in all Figs. 5, 6, and 7. 
Contour is not drawn. The Thomas-Fermi density is l/27rft 2 = 63.7. The energy is E = 0.814. 

FIG. 6. Same as Fig. 5 for E = 0.794. 

FIG. 7. Same as Fig. 5 for E = 0.746. 

FIG. 8. The exact scar strength A (full line) and the semi-classical one calculated with 5 orbits 
(dashed), as a function of energy, at three points of coordinate space. Bottom: x = 0, y = —0.65, 
emphasizing orbit no. 1. Middle: x = 0.8, y = 0, emphasizing orbit no. 2. Top: x = 0.75, y = 0.75, 
emphasizing orbit no. 3. 
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